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Abstract 

The chiral phase boundary of strong matter is determined in the T— \i plane from the chiral quark 
model, applying a non-perturbatively renormalised treatment, involving chains of pion-bubbles 
and 1-loop fermion contributions. In the absence of explicit symmetry breaking the second order 
portion of the phase boundary and the location of the tricritical point (TCP) are determined 
analytically. Sensitivity of the results to the renormalisation scale is carefully investigated. The 
softening of the cr-pole near the second order transitions is confirmed. 
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1 Introduction 



The chiral phase structure of QCD at nonzero quark density receives renewed interest 
with the advent of new lattice techniques [1,2,3,4,5,6]. Recent investigations concentrate 
on the transformation of the temperature driven crossover into first order transitions with 
increasing bario-chemical potential. In the present letter we shall address the same issue. 

The lattice approach is an exact tool, but requires thorough investigation of finite quark 
mass and lattice spacing effects. In particular, the chiral limit with zero mass pions is 
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beyond the access of present numerical QCD techniques (see, however [7]). An interesting 
attempt for a "phenomenological" interpretation of lattice data in the hadronic phase was 
recently proposed. A hadronic resonance model with bag-motivated quark mass depen- 
dence of the hadron spectra describes fairly well the thermodynamics of lattice QCD with 
unnaturally heavy pseudoscalar mesons in the low temperature phase [8,9] (see also [10]). 

Extensive literature exists on the phase structure of effective theories matched to phe- 
nomenological data. Following the pioneering work of Barducci et al. [11] detailed stud- 
ies were published in Nambu-Jona-Lasinio models [12,13,14] and in the linear a model 
[15,16,13]. 

These investigations do not go beyond some version of the quasi-particle approach, e.g. the 
thermodynamical potential is approximated by a sum of free particle contributions with 
masses obtained from some sort of selfconsistent gap equations. The renormalisation of 
the selfconsistent set of equations of state (EoS) and propagators represents a notoriously 
difficult problem. In particular in the broken phase the UV divergences of the selfenergies 
have temperature and /i-dependent coefficients, also depending on the chiral condensate. 
In general, subtraction of such terms defines oversubtracted renormalisation schemes. In 
higher orders of the perturbative solution one expects that resummed or improved expan- 
sions demonstrate the cancellation of the apparent temperature//z/condensate-dependent 
counterterms [17]. Despite promising recent progress [18,19] no general enough implemen- 
tation is available to date for obtaining a renormalised set of equations in which resummed 
propagators could be used. Therefore it is hard to assess the accuracy of the predictions 
obtained from quasi-particle descriptions. 

In the present letter we extend our recent treatment of the pure 0(N = 4) model in the 
approximation of large number of Goldstone bosons (N) [20] to the chiral quark model 
at finite /z. The effect of the constituent quarks will be represented in this letter by 1- 
loop contributions. We shall find fermion effects which are of 0{l/\/~N). A resummed 
renormalisation scheme will be proposed for the coupled set of EoS, and the n and a 
propagators. Finally, the renormalised equations will be analysed in the chiral limit, with 
Goldstone's theorem being consistently obeyed. Since N = 4 is only moderately large, the 
fermionic contributions will not be considered in the numerical solution as perturbations 
of the leading large N solution. They rather will be included fully into the equations, and 
we shall search for a new (non-iterative) solution. The phase boundary will be fully traced 
in the T — \i plane. 

The parameters of the model will be fixed with reference to the T = mass and width 
of the a particle [21]. It turns out that the renormalisation scale influences very strongly 
the shape of the spectral function in the sigma channel. The reason for this phenomenon 
is found in the appearance of extra unphysical poles on the physical Riemann sheet in a 
certain range of the renormalisation scale. After restricting the investigation to the phys- 
ically allowed region, we decided to scan through a large domain of the renormalisation 
scale with the aim to see the sensitivity of the phase boundary in the T — //-plane to its 
choice. We find that T c (/i = 0) is fairly close to the range of critical temperatures mea- 
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sured in two-flavor QCD with different fermion discretisations. The coordinates of TCP 
are in good qualitative agreement with its earlier determinations from effective models 
[13], implying T TC p < 70 MeV independently of the renormalisation scale. 



2 The model and its equation of state 



In the definition of the model we display explicitly the flavor factors: 
L [a, 7r a , V] = L M [<t, tt"] + L F [a, 7r a , ^] + 5L ct [a, ir a , V] , 

where 
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SL ct [a, 7r a , ip] is the countertem Lagrangian. In the large N limit one has Nf — y/N. The 
above form refers explicitly to the broken symmetry phase with the condensate giving rise 
also to the mass of the constituent fermion, which stays finite in the large N limit: 



a(x) -> VN$ + a(x), 



m„ 



(3) 



It is noteworthy that the linear a model has a second independent quartic coupling for 
Nf > 2, which we fix for simplicity at zero, since we wish to evaluate the expressions 
derived in the large N approximation for N = 4. Also, in this letter we work in the chiral 
limit, that is no explicit symmetry breaking term is introduced into the Lagrangian. The 
chemical potential for fermions is introduced into the fermionic piece of the Lagrangian 
by the replacement: d t — > d t — ifJ,. 



(4) 



The form of EoS to 0(1/ \/ N) with A" c colored quarks is the following: 

9 N c,7 



N<5> 



m2 + i* 2 + ^< rt ">" 



+ 



(W) = 0. 



It is obvious that the fermionic contribution is 0(1/ v N). The fermion tadpole will be eval- 
uated using the tree-level mass m q . The subscript M at the expectation value of the com- 
posite operator n 2 (x) refers to an effective pion tadpole computed with propagator mass 
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M, which sums the contribution of all super-daisy type diagrams [22] . The corresponding 
pion propagator is determined from the Schwinger- Dyson equation as iG~ l {p 2 = M 2 ) = 0: 

iG~\p) = p 2 -m 2 - V - ^KOm - Bi*(p). (5) 

Here is the contribution of the fermionic bubble to the self-energy of the pion. By 
an algebraic transformation of its integrand it can be related to the fermionic tadpole as 
follows: 

where I F {p,m q ) is the scalar bubble integral with external momentum p, but evaluated 
with the fermion mass and with Fermi-Dirac distribution, f F (uj) = l/(exp(/3(u; =F/i)) + 1): 

I F (p, m q ) = Ip{p, m q ) - I F {p, m q ), (7) 
I° F (p,m q ) = J (27r)4(A;2 _ m 2 + ie ^ k _ p) 2_ m 2 + ie y ( 8 ) 



I F (P, m q) 



d 3 k 1 

(27r) 3 4wi^ 2 



Po - ooi - uj 2 + ie po + uo\ + uo 2 + 
p - uoi + lu 2 + po + - cj 2 + 



(9) 



Substituting (6) into (5) and making use of (4) one finds that for p 2 = the inverse 
propagator vanishes. This means that Goldstone's theorem is fulfilled and a zero mass 
pion can be used consistently for the evaluation of the pion tadpole in (4). 

Next, we evaluate EoS, for instance, with dimensional regularisation. After MS sub- 
traction, which corresponds to a perturbative renormalisation of A with a counterterm 
proportional to g 4 , one arrives for EoS at 

m 2 + \ ($ 2 + T B (0)) - ^T F (m q ) = 0, (10) 

where T a (x), a = B,F stands for the bosonic/fermionic tadpole contribution with mass 
x: 



T a (x) = —z In -j^- - T;(x), T^x) = — ^ d W V^/„H, 

/ F (a;)^[/ F »+/ F »] ) / fl ( fa; ) = _ s J_. (11) 

Here M 0a are the renormalisation scales, which can be chosen different for the fermionic 
and the bosonic tadpoles. Below we shall characterize the relationship of the two scales by 
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r\ = 1ii(Mob/Mof), which introduces an extra finite subtraction into the renormalisation 
scheme. We shall check if substantially different phenomenology can be achieved when 
this freedom is exploited. 



3 Analytical determination of the TCP 



When looking for the location of continuous phase transitions in the T — [i plane, one can 
study the left hand side of (10) for small field values and in power series of the fugacity, 
exp(/i/T). One expands the distribution function for fermions and antifermions in power 
series, and changes the order of the integration and the summation, doing the integration 
first. One expands the resulting modified Bessel functions up to linear order in m q and 
eventually performs the remaining sum. Remarkably, the logarithmic dependence on $ is 
cancelled, and a power series in <3> is the result. The sign for a continuous phase transition, 
the vanishing of the constant term is given (after substituting N = 4) by the equation 

m 2 + A T 2 _ gi Nc (Li 2 (-e^) + Li 2 (-e-^)) = 0. (12) 
In particular for ji = one finds a shift downwards in T c due to the fermionic fluctuations: 

^) Tl = 0. (13) 

For small fi the coefficient of the quadratic part is positive, therefore T c is consistently 
interpreted as a critical point. The second order line ends at T = T TC p, \i = Htcp, where 
also the quadratic coefficient changes sign: 




A g 4 N c 
6 4tt 2 
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n=0 



cT TC p 

111 w 



0, (14) 



(ln(c/2) = 1 — Je + v)- The first order line for larger \i values was determined numerically 
directly from (10). Also the perfect agreement of the semi-analytically determined second 
order portion of the phase boundary with the exact one was verified numerically. 

The parameters A, g, $, M 0a , m 2 should be determined from the T = data. One 
requires the EoS to give $ = f n /2, with f w = 93 MeV. Then the value of m q (T — 0) = m q0 
determines g. We have chosen m g0 = m N /3 ~ 312.67, so that g = 6.72. The EoS at T = 
provides one relation among m 2 , A and the renormalisation scales M 0a , from which m 
can be determined once the other three parameters are known. One would like to fully 
fix these three parameters M 0a and A, for instance, by studying the spectroscopy of the a 
particle [20]. The relationship between the pole and the spectral function will be discussed 
in section 4. Here we elaborate on the features of the phase boundary when scanning 
through a large domain of M 0B and i] for fixed value of A = 400. 

In Fig.l (l.h.s.) the dependence of the critical temperature on M b is shown at \i = for 



5 




Fig. 1. L.h.s.: the variation of T c (fi = 0) (upper curves) and of Ttcp (lower curves) with Mqb for 
r] = and —0.25 (for the definitions see the text). The arrows on the abscissa separate the region 
of stable (right) and unstable (left) solutions for different values of r/. R.h.s.: the dependence of 
T c (/x = 0),Ttcp an d Htcp ° n the constituent quark mass m q Q. 

i] = 0, —0.25. It should be emphasized that by the arguments to be presented in Section 
4 the allowed range of Mqb lies to the right from the arrow, appearing on the horizontal 
axis of the figure. It is clear that in this region one can arrange very easily T c (/x = 0) to 
occur in the region 150-170 MeV, where it is indicated by lattice studies of the two flavor 
QCD with dynamical quarks [23,24]. The temperature of the tricritical point changes 
very mildly and stays for all choices below 70 MeV. Another informative representation 
of the phase structure is to display the dependence of T c (/i = 0),T TC p and Htcp on the 
constituent quark mass, which appears on the r.h.s. of Fig. 1. 

140 f .-i, 1 1 1 1 1 1 1 1 

120 . >.-400 
100 
> 80 

<D 

I— 60 
40 

--- 2 nd order line 
20 — 1 st order line 

spinodal 

1 1 — ^ 1 ' 

100 200 300 400 500 600 700 800 900 1000 

ii B [MeV] 

Fig. 2. The T — \i phase diagram in the chiral limit for Mqb = 886 MeV and rj = 0. 

We draw the complete phase diagram in Fig. 2 with m q0 = m N /3 and rj = for M 0B = 886 
MeV. The left spinodal is simply the continuation of the second order phase boundary, the 
right spinodal is found by locating the turning point of the two-valued order parameter 
curve as a function of \x for fixed T. We notice that the metastable region between the 
spinodal curves becomes narrower with increasing M 0B - In between we draw the phase 
boundary by finding the degeneracy point of the thermodynamical potential. When com- 
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pared with the features of the previous (non-renormalised) investigations in effective mod- 
els [13], one realizes that the theoretical consistency achieved in the present treatment 
does not produce qualitative differences. If we tune rj in such a way that T c (/i = 0) in- 
creases then TCP moves further down and rightwards on the phase diagram. Our phase 
boundary drops much steeper than the present lattice estimates [6] performed with rather 
heavy quarks and is quite close to the freeze-out curve of Cleymans and Redlich [25]. 



4 Poles in the a channel and their relation to the spectral function 



In the present letter we need the spectral function in the a channel only for T = 0, 
that is for discussing the sensitivity of its pole structure to the selection of the couplings, 
especially the renormalisation scales of the effective model. A detailed analysis of the 
features of the finite temperature pole trajectory and the phenomenon of the threshold 
enhancement is postponed to a forthcoming publication. The large N form of the inverse 
propagator of the a particle [20] is now modified by the contribution of the fermionic loop: 



iG~\p) =P 2 - m 2 - ^ - ^(^ a )M - S chain ( P ) 

- ^^(V'V') + 2 ^^ c (4m 2 - p 2 )I F (p, m q0 ) + counterterms. (15) 

Here S c h a in{p) represents the sum of the infinite sequence of pion bubble contributions to 
the a self energy. This equation simplifies if (4) is taken into account to 



A 2g N c 
iG~ 1 (p) = p 2 - -$ 2 - S chain {p) H j=^{Am 2 - p 2 )I F {p, m q0 ) + counterterms. (16) 



Fermionic counterterms are needed for the perturbative renormalisation of the coefficients 
of p 2 and of <3> 2 , both receiving divergent contribution from If(p, m q o)- Using MS sub- 
traction the fermionic part of the coupling constant counterterm which is 0(g 4 ), proves 
to be the same as the one encountered in the renormalisation of the EoS. Additional 
bosonic counterterms form a power series in A, and are needed to renormalise the chain 
of bubbles S c h a in- This series can be formally resummed, resulting in a non-perturbative 
renormalisation of the meson sector [20] . The relation between the bare and renormalised 
coupling constant is: 

Abarc = - ; i + ^-T= b div, b div = — (-- + 7 - 2 - \h(Att)) . (17) 



Summing up the renormalised terms of the bubble chain, the final renormalised a propa- 
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gator at rest (p = 0) reads as follows: 



tG- 1 ( Po ,0)=p 2 - ^ 2 (l - ^I b ,r(po,0) 



-i 



+ 



2g 2 N c 



N 



(4mJ - pI)I f ,r(po, m q0 ). (18) 



The spectral function is obtained by putting p — > p + ie, and taking the imaginary part 
along the physical real axis: 



Pa(po) = -Im [iG a (p Q + ie, 0)] . 

7T 



(19) 



Below we explore how the a pole structure in the lower halfplane is reflected in the shape 
of the spectral function for a broad range of the values of the renormalisation scale M$ B 
as it appears in Fig. I. An analytic continuation of (18) is understood via the T = 
expressions of I a ,n, which for complex values of po can be written explicitly as 



JjU(Po»"Ho) = 



16tT 2 
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90 



(20) 



The above formulae correspond to analytical continuations onto specific Riemann-sheets 
characterised by two integers (fci,/^). For example k\ — 1, k 2 = 1 and fci = 1,&2 = — 1 
corresponds to an analytical continuation across the real axis above and below the two- 
fermion threshold, respectively. We found roots of the right hand side of (18) for all 
fci > 0, k 2 < integers. For a wide range of A values the real parts of all poles have nearly 
the same value and lie slightly below the two-quark threshold. Therefore one expects the 
spectral function to have the maximum of its "resonant" part at po tma x ~ Re(j9 ,po«e)- 

This expectation turned out to be incorrect below certain values of the renormalisation 
scale, depending on rj. What spoils this feature is the Landau pole structure on the first 
Riemann sheet, which is an indication of the effective nature of the theory. This structure 
sensitively depends on M 0B . Apart from the well known large scale Landau pole(s) on 
the imaginary axis, which decreases exponentially with increasing values of A, new small 
scale poles might emerge. 

The mechanism for their occurrence is the cleanest to illustrate for A = 0. In this case the 
imaginary part of a physical a pole necessarily vanishes. For M QB > M 0L the pole is real 
but its location shifts towards the origin as M 0B decreases. At M 0B = M 0L it reaches the 
origin and collides with its mirror arriving along the real axis from the direction p < 0. 
For M 0B < M l a conjugate pair of poles appears along the imaginary axis. One of them 
signals the instability of the system at very low energy scales, completely invalidating 
our solution in this region of the renormalisation scales. This dipole induced instability 
prevails also for A = 400, restricting the choice of the renormalisation scale to M 0B > M 0L . 
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Fig. 3. The T = spectral function for different allowed choices of Mq B and r] in the chiral limit. 
The variation is restricted to the stability region of the system. 



One finds from (20) Mql = m q oexp(j] + 1) which, with our present numerical values of 
f„ and m q0 , gives M 0L = 849.92 MeV at rj = and M 0L = 661.92 MeV at rj = -0.25. In 
the allowed region M 0B < M 0L , \M Landau \/\p 0tPO i e \ stays always well above unity. There, 
the real part of the pole position and the maximum of the spectral function coincide as 
expected. 



In Fig. 3 we display the emerging spectral functions for M 0B = 850, 886 MeV and rj = 0. 
In addition also the curve for M 0B = 680 MeV and rj = —0.25 is shown. The effect of the 
fermions brings the location of the "cr-peak" closer to the phenomenologically expected 
region: at rj = we find po, P eak < 430 MeV, while the corresponding poles are located 
at Repose — 450 MeV. The Y a /M a ratios still turn out to be lower than expected by 
40-50%. Finally, one might remark that the stability of the broken symmetry phase also 
imposes a lower limit on M 0B , which is, however, lower in our case than the bound, 
described above. 



Two short notes are of order in connection with the thermal evolution of the poles. First, 
near T c all "cr-poles" flow to the origin of the po plane asymptotically along the same curve. 
The real part of this trajectory scales with $(T), the imaginary part with $ 2 (T). This 
means that the inclusion of the fermions qualitatively changes the dynamics of the system, 
since in the pure O(N) model the trajectory reached the origin along the imaginary axis. 
Second, for the scaling exponent of the order parameter, $(T) ~ (T — T c ) /3 ', on the full 
second order line we find its mean field estimate 1/2. For the single point of TCP /3 — 1/4 
was found within the numerical accuracy of its determination, in full agreement with 
the standard result of statistical physics [26]. At TCP only logarithmic corrections are 
expected to this value. 
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5 Conclusions and discussion 



In this letter we have presented a theoretically fully consistent derivation of the phase 
structure in the T — [l plane of strong matter from the chiral a model in the chiral limit 
using an approximate solution valid to 0(1/ VN) when the number of the Goldstone 
bosons (N) goes to infinity A novel feature arising from including the fermions is the 
strong influence of the renormalisation scheme on the shape of the spectral function. It 
is essential to check also the dependence of the excitation spectra on the renormalisation 
scale, which allows to restrict considerably the range for the acceptable values of M 0B . 
The predictions on the phase boundary in the T — \x plane depend on the renormalisation 
scale(s) only smoothly, leading to a unique semi-quantitative characterisation of the phase 
structure. In the allowed range of the parameters one finds T c (fi = 0) in reasonable 
agreement with the results of MC simulations, while one encounters a genuinely low value 
for T TCP ^ 60 - 70 MeV. 

Our approach can be further improved by the dynamical determination of the fermion 
mass from a one-loop Schwinger- Dyson equation. For a realistic comparison with hot 
heavy-ion systems it is important also to include the effect of explicit chiral symmetry 
breaking. Possible dependence on the renormalisation scheme can be optimised by re- 
quiring the best phenomenological description of the scalar-isoscalar spectral function 
for T = 0. These extensions require a careful implementation of the proposed (non- 
perturbative) renormalisation procedure, which will prove more complex than in the 
present case. Our final goal is the application of the large N method to the three-flavor 
SU(3)l x SU(3)r case with realistic constituent quark masses. 
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